功能梯度压电材料圆环位移和电势边值问题的精确解

圆环形的压电材料器件在智能结构中得到了广泛的应用。本文推导了横观各向同性功能梯度压电材料圆环在内、外边界上给定位移和电势情况下的一般解。极化方向在圆环的半径方向，材料常数的梯度方向也设定在半径方向，并可表示为半径r的幂，本构关系为线性。然后推导了压电圆环外壁固定、接地，内壁沿垂向有一微小位移、电势分别为余弦分布和均匀分布的问题的精确解，并计算了该问题在这两种电势情况下产生的无量纲形式的径向和环向位移、电势、应力及电位移沿径向分布的数值结果。计算中考虑了不同的材料梯度，以及内壁的位移与电势的不同比例。     

功能梯度材料板件三维分析的半解析梯度有限元法

将半解析有限元与梯度有限元相结合,形成一种半解析梯度有限元来求解功能梯度材料板件问题。该方法兼有有限元法的适应性强、程序统一,半解析有限元法的节省单元与计算工作量,梯度有限元法的适应构件内部材料性能任意梯度分布等特点,并实现用一维数值计算给出构件三维分析结果。算例分析表明了方法的精度、功能与上述特点,充分揭示了功能梯度材料板件力学响应的三维形态。半解析梯度有限元法可推广应用到其他功能梯度材料面结构的各类分析中。     

The transient piezothermoelastic problem of a thick functionally graded thermopiezoelectric strip due to nonuniform heat supply

This paper is concerned with the theoretical treatment of the transient piezothermoelastic problem involving a thick functionally graded thermopiezoelectric strip due to nonuniform heat supply in the width direction. The thermal, thermoelastic and piezoelectric constants of the strip are assumed to vary exponentially in the thickness direction. The transient two-dimensional temperature is analyzed by the methods of Laplace and finite sine transformations. We obtain the exact solution for a simply supported strip under the state of plane strain. Some numerical results for the temperature change, the displacement, the stress and electric potential distributions are presented in figures and table. Furthermore, the influence of the nonhomogeneity of the material and that of the electric boundary conditions are investigated.     

功能梯度压电材料反平面裂纹问题

基于三维弹性理论和压电理论，导出了材料系数在横观各向同性平面内梯度分布的压电体的状态方程，进而对材料系数指数函数规律分布的半无限大压电体中的反平面裂纹问题进行了求解，利用Fourier变换给出了半无限大压电体中位移，应力，电势及电位移的解析表达式，并求得了裂纹尖端的应力强度因子和电位移强度因子，分析了不同的非均匀材料系数及几何尺寸对它们的影响。     

Dynamic stress of a circular cavity buried in a semi-infinite functionally graded piezoelectric material subjected to shear waves

The paper presents a theoretical method to investigate the multiple scattering of shear waves and dynamic stress around a circular cavity in a semi-infinite functionally graded piezoelectric material. The analytical solutions of wave fields are expressed by employing wave function expansion method and the expanded mode coefficients are determined by satisfying the boundary conditions of the cavity. Image method is used to satisfy the free boundary condition of the semi-infinite structure. According to the analytical expression of this problem, the numerical solutions of the dynamic stress concentration factor around the cavity are presented. The effects of the piezoelectric property, the buried depth of the cavity, the incident wave number and the nonhomogeneous parameter of materials on the dynamic stress around the cavity are analyzed. Analyses show that the piezoelectric property has great effect on the dynamic stress in the region of intermediate frequency and the effect increases with increasing wave number. When the nonhomogeneous parameter of materials is less than zero, it has less influence on the maximum dynamic stress around the cavity; however, it has greater influence on the distribution of the dynamic stress around the cavity. When the nonhomogeneous parameter of materials is greater than zero, it has greater influence on both the maximum dynamic stress and the distribution of dynamic stress around the cavity, especially in the case that the buried depth is comparatively small.     

Exact solutions for radial deformations of a functionally graded isotropic and incompressible second-order elastic cylinder

We find closed-form solutions for axisymmetric plane strain deformations of a functionally graded circular cylinder comprised of an isotropic and incompressible second-order elastic material with moduli varying only in the radial direction. Cylinder's inner and outer surfaces are loaded by hydrostatic pressures. These solutions are specialized to cases where only one of the two surfaces is loaded. It is found that for a linear through-the-thickness variation of the elastic moduli, the hoop stress for the first-order solution (or in a cylinder comprised of a linear elastic material) is a constant but that for the second-order solution varies through the thickness. The radial displacement, the radial stress and the hoop stress do not depend upon the second-order elastic constant but the hydrostatic pressure and hence the axial stress depends upon it. When the two elastic moduli vary as the radius raised to the power two or four, the radial and the hoop stresses in an infinite space with a pressurized cylindrical cavity equal the pressure in the cavity. For an affine variation of the elastic moduli, the hoop stress in an internally loaded cylinder made of a linear elastic isotropic and incompressible material at the point is the same as that in a homogeneous cylinder. Here R_in and R_ou equal, respectively, the inner and the outer radius of the undeformed cylinder and R the radial coordinate of a point in the unstressed reference configuration.     

The Stress Response of Functionally Graded Isotropic Linearly Elastic Rotating Disks

The purpose of this research is to investigate the effects of material inhomogeneity on the response of linearly elastic isotropic solid circular disks or cylinders, rotating at constant angular velocity about a central axis. The work is motivated by the recent research activity on functionally graded materials (FGMs), i.e., materials with spatially varying properties tailored to satisfy particular engineering applications. The analog of the classic problem for a homogeneous isotropic rotating solid disk or cylinder is considered. The special case of a body with Young"s modulus depending on the radial coordinate only, and with constant Poisson"s ratio, is examined. For the case when the Young"s modulus has a power-law dependence on the radial coordinate, explicit exact solutions are obtained. It is shown that the stress response of the inhomogeneous disk (or cylinder) is significantly different from that of the homogeneous body. For example, the maximum radial and hoop stresses do not, in general, occur at the center as in the case for the homogeneous material. Furthermore, for the case where the Young"s modulus increases with radial distance from the center, it is shown that radially symmetric solutions exist provided the rate of growth of the Young"s modulus is, at most, cubic in the radial variable. It is also shown for the general inhomogeneous isotropic case how the material inhomogeneity may be tailored so that the radial and hoop stress are identical throughout the disk. This revised version was published online in August 2006 with corrections to the Cover Date.     

含压电层的功能梯度柱壳的自由振动分析

基于三维弹性理论，导出了带有压电层的圆柱形梯度壳的动力学方程以及相应的边界条件，用幂级数展开法得到了求解该圆柱形梯度壳自由振动的三维精确公式．通过实例模型求解了该壳体的自由振动的固有频率；分析了不同电学边界条件对壳体的振动频率的影响。结果可评估各种近似理论解和数值解的正确性。     

Effect of material nonhomogeneities on the HRR dominance

This work studies the asymptotic stress and displacement fields near the tip of a stationary crack in an elastic–plastic nonhomogeneous material with the emphasis on the effect of material nonhomogeneities on the dominance of the crack tip field. While the HRR singular field still prevails near the crack tip if the material properties are continuous and piecewise continuously differentiable, a simple asymptotic analysis shows that the size of the HRR dominance zone decreases with increasing magnitude of material property gradients. The HRR field dominates at points that satisfy |α⁻¹ ∂α/∂x_δ|1/r, |α⁻¹ ∂²α/(∂x_δ ∂x_γ)|1/r², |n⁻¹ ∂n/∂x_δ|1/[r|ln(r/A)|] and |n⁻¹ ∂²n/(∂x_δ ∂x_γ)|1/[r²|ln(r/A)|], in addition to other general requirements for asymptotic solutions, where α is a material property in the Ramberg–Osgood model, n is the strain hardening exponent, r is the distance from the crack tip, x_δ are Cartesian coordinates, and A is a length parameter. For linear hardening materials, the crack tip field dominates at points that satisfy |E_tan⁻¹ ∂E_tan/∂x_δ|1/r, |E_tan⁻¹ ∂²E_tan/(∂x_δ ∂x_γ)|1/r², |E⁻¹ ∂E/∂x_δ|1/r, and |E⁻¹ ∂²E/(∂x_δ ∂x_γ)|1/r², where E_tan is the tangent modulus and E is Young’s modulus.     

End Effects in Anti-plane Shear for an Inhomogeneous Isotropic Linearly Elastic Semi-infinite Strip

The purpose of this research is to further investigate the effects of material inhomogeneity on the decay of Saint-Venant end effects in linear isotropic elasticity. This is carried out within the context of anti-plane shear deformations of an inhomogeneous isotropic elastic solid. The mathematical issues involve the effects of spatial inhomogeneity on the decay rates of solutions to Dirichlet or Neumann boundary-value problems for a second-order linear elliptic partial differential equation with variable coefficients on a semi-infinite strip. In previous work [1], the elastic coefficients were assumed to be smooth functions of the transverse coordinate so that the material was inhomogeneous in the lateral direction only. Here we develop a new technique, based on a change of variable, to study generally inhomogeneous isotropic materials. The governing partial differential equation is transformed to a Helmholtz equation with a variable coefficient, which facilitates analysis of the influence of material inhomogeneity on the diffusion of end effects. For certain classes of inhomogeneous materials, an explicit optimal decay estimate is established. The results of this paper are applicable to continuously inhomogeneous materials and, in particular, to functionally graded materials. This revised version was published online in August 2006 with corrections to the Cover Date.